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Abstract 



The mSUGRA model predicts the polarization of the tau coming from the stau to bino decay 
in the co- annihilation region to be +1. This can be exploited to extract this soft tau signal 
at LHC and also to measure the tiny mass differences between the stau and the bino LSP. 
Moreover this strategy will be applicable for a wider class of bino LSP models, where the 
lighter stau has a right component at least of similar size as the left. 



1. Introduction 

The mininial supersyinmetric standard model (MSSM) has been the most popular extension 
of the standard model(SM) for three reasons. It provides a natural solution to the hierarchy 
problem of the SM and a natural candidate for the cold dark matter in terms of the lightest 
superparticle(LSP) along with the unification of gauge couplings at the GUT scale. In 
particular there is a great deal of interest in the minimal supergravity (mSUGRA) model 
as a simple and well-motivated parametrization of the MSSM. This is described by the four 
and half parameters [1] 

"^i/2,'"^o, ^0, tan/? and sgn(/i), (1) 

the first three representing the common gaugino and scalar masses and trilinear coupling at 
the GUT scale. The tan/3 stands for the ratio of the two Higgs vacuum expectation values, 
while the last one denotes the sign of the mixing parameter n between them. The magnitude 
of fi is fixed by the radiative electroweak symmetry breaking condition. 

Astrophysical constraints on dark matter(DM) require the LSP to be colorless and neu- 
tral, while direct DM search experiments strongly disfavor sneutrino LSP. Thus the favored 
candidate for LSP in the MSSM is the lightest neutralino, 

X? = NnB + NuWs + iVis^i + iVi4^2. (2) 

In the mSUGRA model the Xi is dominantly bino(-B) over most of the allowed parameter 
space. Since bino has no gauge charge, it can pair-annihilate mainly via sfermion exchange. 
The large sfermion mass limits from LEP [2] makes this annihilation process inefficient, 
leading to an overabundance of DM over most of the mSUGRA parameter space. There are 
essentially two narrow strips of cosmologically compatible DM relic density [3] at the opposite 
edges of the parameter space, corresponding to the large and the small mo boundaries, called 
the focus point and the stau co-annihilation regions respectively [4, 5]. There is also a narrow 
strip in the middle called resonant annihilation region, but only at very large tan/3(~40). 
Out of these only the stau co-annihilation region is compatible with the muon anomalous 
magnetic moment constraint [6], and also the only one which can be completely covered 
at the LHC. Therefore the stau co-annihilation region is a region of special interest to the 
SUSY search programme at the LHC. In particular one is looking for a distinctive signature 
which will identify the SUSY signal at the LHC to this region and also enable us to measure 
the tiny mass difference AM between the co-annihilating superparticles, which is predicted 
to be ~5% by the DM relic density constraint. A distinctive feature of this region is that a 
large part of the SUSY cascade decay occurs via 

n ^ rxl (3) 

leading to a r lepton along with the canonical missing Et{I^j^)- Unfortunately, the r leptons 
are very soft because of the small mass difference between fi and Xi- But fortunately the 



r polarization is predicted to be very close to +1 in the mSUGRA model. In this work we 
hope to show with the help of a generator level Monte Carlo simulation that the positive 
polarization(PT- = +1) of this r signal can be exploited to extract it from the negatively 
polarized r (Pr = — 1) background as well as the fake r background from hadronic jets. 
Moreover, the steep pt dependence of the soft r signal will provide a distinctive signature 
for the co-annihilation region as well as a measure of the tiny mass difference between the 
co-annihilating particles. 

In section 2 we briefly describe how to use the r polarization. In section 3 we describe 
the stau co-annihilation region with the help of an illustrative point in the parameter space 
of this region. We also identify the main channels of SUSY cascade decay. In section 4 we 
discuss the result of the simulated SUSY signal for this point. We show how the polarization 
cut retains most of the Pr = +1 signal r-jets while suppressing Pr = —1 background r-jets 
and practically eliminating the fake r jets. We also estimate the significance level of the 
signal and the AM measurement from the slope of the soft r jet signal. We conclude with 
a summary of our results in section 5. 

2. Using T Polarization 

The best channel for r identification is its 1-prong hadronic decay channel, accounting for 
50% of its decay width. Over 90% of this comes from 



n' 



=z/(12.5%), p±z/(26%), afu{7.5%), (4) 



where the branching fractions for ir and p include the small K and K* contributions respec- 
tively, which have identical polarization effects [2]. The CM angular distribution of r decay 
into vr or a vector meson v{p, ai) is simply given in terms of its polarization as, 

Y^l a = -{I + Pr cosO) 

L TT a cos u 2 

1 dV^L^T lm1,ml 

Y^i — a = -^—TT^y^^^rCose) (5) 

L V a cos u m^ + Zmf, 

where L,T denote the longitudinal and transverse polarization states of the vector meson. 
The fraction x of the r laboratory momentum carried by its decay meson i.e. the (visible) 
r-jet, is related to the angle 6 via 

1 771 

x = -(1 + cos^) + ^(1 - cos^), (6) 

in the collinear approximation(p^ ^ m^). It is clear from eqs. 5 and 6 that the relatively 
hard part of the signal {Pr = +1) T-jet comes from the vr, p^ and aiL contributions; while 
for the background(P^ = ~1) T-jet it comes from p^ and aix contributions [7]. This is the 
important part that would pass the pt threshold for r jets. 



Now the pt and aix decays favor even sharing of the momentum among the decay pions, 
while the pi and aiL decays favor uneven distributions where the charged pion carries either 
very httle or most of the momentum. Thus plotted as a function of the momentum fraction 
carried by the charged pion 

i?=^^, (7) 

Vr—jet 

the longitudinal p and ai contributions peak at very low or very high R(~ 0.2 or ~ 0.8), while 
the transverse contributions peak in the middle [7, 8]. The low R peak of the longitudinal 
p and ai contributions are not detectable because of the minimum pt requirement on the 
charged track for r identification {R >0.2). Now moving the R cut from 0.2 to 0.8 one cuts 
out the transverse p and ai peaks while retaining the detectable longitudinal peak along with 
the single tt^ contribution. Thanks to the complementarity of these two sets of contributions, 
one can effectively suppress the former while retaining most of the latter by a simple cut on 
the ratio 

R > 0.8. (8) 

Thus one can suppress the hard part of the r-jet background(PT- = —1), while retaining most 
of it for the signal {Pr = +1), even without separating the different meson contributions from 
one another [8] . This is a simple but very powerful result particularly for the hadron colliders, 
where one cannot isolate the different meson contributions to the r-jet in eq. (4). This has 
been used to enhance the Pr = +1 signal of charged Higgs boson at Tevatron and the LHC 
[8, 9]. It can be also used with effect in the investigation of the Pr = +1 SUSY signal coming 
from eq.(3) at Tevatron and the LHC [10] as well as the ILC [11]. 

This is the first application of the r polarization effect, however, in the context of probing 
the stau co-annihilation region. As we shall see below the simple polarization cut of eq. (8) 
helps to suppress not only the {Pr = —1) T-jet backgrounds, but simultaneously the fake 
r jet background from hadronic jets as well. We shall also see that the conclusion of this 
section remains true even after we add the non-resonant contributions to the r decay of 
eq.(4). 

3. Stau co-annihilation region of the mSUGRA model 

For illustration we have chosen one point in the stau co-annihilation region of the mSUGRA 
model. We expect the results to hold equally for any other point in this region. The input 
mSUGRA parameters are shown in Table-1 along with the resulting weak scale superparticle 
spectrum calculated using ISAJET(v7.74)[12]. For simplicity we have taken Aq = and 
positive sign for p since our results are not sensitive to them. The corresponding LHC cross 
sections for the three strong processes at LO are, 

^i~99, M, m + qq*) = (0.46, 2.4, 1.3)pb. (9) 



These are evaluated with the CTEQ3L[13] structure functions setting a common factorization 
and renormahzation scale at the average mass of the superparticle pair. But the results are 
insensitive to either of these choices. We have also checked that changing from CTEQ3L to 
the more recent structure functions of CTEQ5L only change the total signal cross-section of 
eq.(9) from 4.16 pb to 4.11 pb, although there are somewhat larger variations for the three 
individual processes. The largest contribution to the third process of eq.(9) comes from qq 
-^ qq because of the dominance of the valence quark flux over gluon. Therefore, we expect 
uncorrelated production of singlet (S) and doublet (D) squarks in all the three processes i.e. 
SS:SD:DD=l/4:l/2:l/4. Thus 3/4 of the events contain one or two doublet squarks. The 
doublet squarks mainly undergo cascade decay via the W dominated chargino and neutralino 
states Xi ? X25 ^^^o, 



xf 



I ~ 

TTi ■ 



r'rxl 



(10) 
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Table 1: Masses of superparticles(in GeV), for an illustrative mSUGRA point in the 

stau-coannihilation region. 



mo 


"^1/2 


tan/9 


9 


qi 


qR 


n 


T-2 


xt 


xt 


X? 


X^ 


A3 


Xa 


223 


400 


40 


942 


~884 


-856 


172.5 


356.8 


307.5 


526 


163 


307 


511.5 


526 



Thus one expects to see one soft r in most cascade decays, sometimes accompanied by a 
relatively hard r' . The lighter stau 



Ti = T~R sin Of + tl cos Of, 



(12) 



is dominated by tr, since there is no SU(2) contribution to the RGE. The resulting polar- 
ization of the soft r from fi -^ rx is given by [10, 14], 



P. 



"11 



On 



^fi)^ 



r(r^) - r(r^) 

r(r^) + r(rz.) (afi)2 

2g 
-i^Nii tan^^y sin6'f — 

v2 

9 



Af 






\^myy COS P 

gm 



iVi3 cos Of 



—7= [Ni2 + Nil tan Ow] cos Of - 

V2 y/Zmw cos p 



iVis sin6'T 



(13) 



The first subscript of ai{ refers to the ri and the second to the neutralino Xi- Thus the 
dominant term is afi ~ — ^A^^n tan^^^sin^f, implying P^ ^ +1. In fact in the mSUGRA 
model there is a cancellation between the sub-dominant terms, so that one gets Pr -0.9 
throughout the allowed parameter space [10]. In the stau co- annihilation region of our 
interest Pr -0.98. So one can safely set Pr =+1. 

The polarization of the hard r' from eq.(ll) is obtained from eq.(13) replacing ail by 
ai2 ■ The dominant contribution comes from 0^2 — ^^22 cos^f , implying Pr' ~-l. There is 



a similar cancellation of the sub-dominant contributions, leading to P,-' <-0.98 in the stau 
co-annihilation region. Thus one can set Pt-/=-1. 

Finally it should be noted that probing the stau co-annihilation region via the r polar- 
ization cut of eq.(8) should be applicable in a wider class of MSSM with a B LSP. Since the 
coupling of B to tr is 2 times larger than to tl in eq. (13) one gets Pr ~ 0.6 as long as the 
fi has a right component as least of similar size as the left component [15]. In that case the 
R>0.8 cut will keep at least half of the soft tau signal. 

4. Event simulation, results and discussion 

The SUSY signal events are generated using the event generator PYTHIA(v6.23) [16], which 
simulates superparticle pair production and cascade decay for the spectrum shown in Table 
1. The generated r leptons are then passed through the TAUOLA^ package [17] to simu- 
late r decay, which includes the effect of r polarization in the hadronic decay channels of 
eq.(4) along with the small non- resonant contribution. The generated events are then passed 
through the CMS JET package [18] for jet reconstruction. The jets are constructed using the 
cone algorithm in CMSJET with a cone size of AR =0.5. The kinematic cuts for jet recon- 
struction are Et >15 GeV and \r]\ <4.5. Finally, the missing ExiJ^rp) is reconstructed by a 
vector summation of calorimetric energies. 

The events are then subject to the following selection cuts for triggering and suppression 
of SM background: 

Number of jets > 2 with E^'-' > 100 GeV, \r]j,J < 4.5 

^j, > 250 GeV. (14) 

From the events passing this cut we select those containing at least one r jet at the generator 
level. We then require the softest r jet in each event to satisfy 

15 GeV < p^-^"* < 40 GeV. (15) 

This is our sample of soft r jet events. We try to simulate the efficiency of tracker isolation 
on the r jet by requiring that (i) it has one and only one charged track (leading track) of 
p^^ >6 GeV within narrow signal cone of ARg = 0.1 measured with respect to its calorimetric 
energy deposit, and (ii) there is no other charged track in a surrounding isolation cone of 
ARj = 0.4 with p^ >3(1) GeV. The p^ >1 GeV isolation cut ensures higher purity of r-jets 
at the cost of a lower efficiency of r identification, as genuine r jets can be often accompanied 
by such a charged track in the environment of the LHC. 

It should be noted here that in a full simulation the tracker isolation cut is supplemented 
with calorimetric cuts for complete r identification [19], which is beyond the scope of the 



^Thc current version of TAUOLA cannot simulate r decay when it comes from sparticle decay. One of 
the authors(MG) has modified it to include these decay modes. 
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present work. Since the efficiency of these supplementary cuts is quite high (~0.8) we shall 
simply assume it to be 1. We know from full simulation studies[19] that r-jets identified 
via tracker isolation and calorimetric cuts match very well (in direction and energy) with the 
generator level r-jets. Since we are unable to incorporate the calorimetric cuts, however, 
we shall work only with generator level r-jets, satisfying the tracker isolation, as in the case 
of ref. [20]. We shall add to these the fake r-jets coming from the generator level hadronic 
jets, which satisfy the tracker isolation cut. Thus we approximately incorporate the tau 
identification efficiency and purity in our analysis. Working with generator level r-jets allows 
us to tag it to its 'mother'. Thus we can separate the soft r jet signal of the stau decay 
of eq. (3) from the background coming from the various other sources in the cascade decay; 
and look at the effect of the polarization cut of eq.(8) on each of them. This will help us to 
understand the major contributors to the soft r jet background. We shall also consider the 
largest SM background to the soft r-jet signal, coming from ti and W+multijet channels. 

Table 2: Number of events for the process qg with and without tracker isolation. "Moth- 
ers" of each r jet, shown below the 3rd row are found by direct tagging in the event generator. 
(I) no tracker isolation (II) with tracker isolation p^ >3 GeV and p!^^ > 6 GeV (III) with 
tracker isolation p^ >1 GeV and p^*" > 6 GeV. Below the 3rd row the parenthetic entry in 
each block is obtained with the R >0.8 cut. 



Selection 


I 


II 


II 


No. of events simulated 


1500k 


1500K 


1500K 


At least one r-jet 
With 15 GeV < p5^-J"* <40 GeV 


431283 
156555 


320207 
108398 


254150 
86137 


h 


87061 
(42259) 


61925 
(35239) 


48415 
(28178) 


xl 


27073 

(10852) 


22813 
(9475) 


19178 
(7913) 


xi 


644 
(345) 


573 

(315) 


509 

(283) 


xl 


712 

(377) 


639 

(353) 


548 
(304) 


w 


22419 
(6770) 


16532 
(5358) 


13033 

(4227) 


z 


787 
(320) 


611 
(265) 


494 
(212) 


h 


2994 
(1189) 


2366 
(1033) 


1925 
(831) 


h 


98 
(40) 


79 

(37) 


59 
(30) 


None of the above 'mothers' 


14767 

(5287) 


2860 
(1265) 


1976 
(939) 



We have generated 1.5 million events for the qg production, which has the largest SUSY 
cross-section of eq.(9) i.e. 2.4 pb. Thus it corresponds to an effective luminosity of about 
600 fb^ . This will ensure that the predicted events have greater statistical accuracy than 
the data, so that our error estimates are primarily controlled by the latter. Table 2 shows the 
effects of the above mentioned cuts on these events. The 2nd row shows that about l/3rd 
of the events contain at least one r-jet, where the suppression factor includes the effects of 
the selection cut (14) and the 1-prong hadronic branching ratio of r lepton. It also shows the 
effect of the tracker isolation cuts in the three columns. We get an efficiency factor of about 
2/3(1/2) corresponding to the p'^ >3(1) GeV on the accompanying charged tracks. The 3rd 
row shows that about 1/3 of these events have the softest r jet in the px range of eq.(15). 
The subsequent rows tag these r jets to their 'mother' and show the effect of the polarization 
cut of eq. (8) in each of these cases in parentheses. We see that about half of these soft r 
jets come from the signal process of eq. (3) having Pr =+1. About 60% of them survive 
the polarization cut. The background is dominated by the Pr =-1 r-jets coming from x% of 
eq. (11), and the W bosons produced in cascade decay. About 30%(40%) of the r-jets from 
W(x2) decay survive the polarization cut (8). The reason for the higher survival probability 
of the r jets from the x\ decay(ll) is that they are quite hard to start with. Hence even the 
relatively soft part of these r-jets, coming from the pi and tt decay channels given in eq. (5), 
survive the modest px cut of eq.(15). The r jets tagged to none of the listed 'mothers' are 
expected to come mainly from b hadron decay. 

Table 3: Estimation of faking efficiency for QCD and SUSY [gg) jet events; QCD jet 
events are for px range of 15-40 GeV, where pt denotes the transeverse momentum in the 
partonic centre of mass frame. 



Selection 


pf >1 GeV, 


p'^' >6 GeV 


pf >3 GeV, 


p'^' >6 GeV 


process^ 


QCD 


SUSY 


QCD 


SUSY 


No. of events simulated 


10*^ 


300K 


10*^ 


300K 


Total no. of jets 


181377 


349989 


827424 


347303 


Total no. of (fake)r jets 

Total no. of (fake)r-jets with R cut 

Faking efficiency 

Faking efficiency with R cut 


6486 

376 

0.034 

0.0015 


9183 
313 
0.026 
0.0009 


123100 
2496 
0.15 
0.003 


38473 

651 
0.11 
0.0019 



Table 3 demonstrates the effect of polarization cut on the fake r jets. It presents the fake 
r jets coming from the hadronic jets produced in SUSY(^^) and QCD processes without 
the selection cuts of eq. (14). It shows that the fake r-jets background from the hadronic 
jets in this SUSY cascade decay are quite large, particularly for the tracker isolation cut of 
p'^ >3 GeV. But it gets practically eliminated by the R cut of eq.(8). The faking efficiency 
of 11(3)% for the isolation cut of pf >3(1) GeV falls to 0.2(0.1)% after the R cut. For 
comparison we also show in this table the faking efficiencies of hadronic jets for 1 million 
simulated QCD events. The faking efficiencies before and after the R cut are seen to be very 
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similar to those of the SUSY events. This is evidently a very powerful result, which is not 
unexpected though. A hadronic jet can fake an one prong r-jet by a rare fluctuation, when 
all but one of the constituent particles (mostly pions) are neutral. Then requiring the single 
charged particle to carry more than 80% of the total jet energy requires a second fluctuation 
which is even rarer. Normally identification of r jets down to pt=15 GeV would be difficult 
at the LHC because of the fake backgrounds, which is alleviated now by the R >0.8 cut. 
Note that the R cut automatically raises the px threshold of the leading track to 12 GeV. 



Table 4: The SM backgrounds from tt and W+multijet processes. In the latter case 
separate simulations are done for the two pT ranges of W, as shown in the table. 



Process— s> 


tt 


W+multijets 






PT =100-300 GeV 


300-lOOOGeV 


No. of events simulated 


3x10^ 


10^ 


10^ 


^y >250 GeV 


15318 


1996 


150073 


Jet cut 


15318 


1996 


150072 


Me// >750 GeV 


9277 


136 


73841 


One r-jet 


1250 


13 


12931 


Mt >50 GeV 


417 


2 


2651 


15 GeV< PT <40 GeV 


150 





602 


From W 


148 





601 


With R cut 


(40) 





(143) 


Raw Cross sections(pb) 


500 


1500 


30 


Number of events for £=10 fb"^ 








Without R cut 


247 


< 1 


180 


With R cut 


(67) 


(<1) 


(43) 



Since the R >0.8 cut practically eliminates the fake r jet background for both the isolation 
cuts, we choose to work with the p^ > 3 GeV cut for two reasons. Firstly it is less demanding 
on the tracker momentum resolution; and secondly it has a higher efficiency for identifying 
genuine r jet, as seen above. We have found that imposing the selection cuts of eq.l4 sup- 
presses the QCD process very effectively in agreement with ref.[20]. Supplementing this with 
the faking efficiency with R>0.8 cut makes the QCD background to the SUSY signal com- 
pletely negligible. The largest SM background comes from the ti and W+multijet processes. 
In order to control these backgrounds we supplement the selection cuts of eqs.(14),(15), with 
two more kinematic cuts, 



M 



eff 



If, 



T 



+ E^i + Ei^ > 750 GeV, 



MT{T-jet,:^^) >50 GeV, 



(16) 



where Mt = J2pr"^'=*^^(l 



COS0(p^ "''^ , 



Ifrp)) and (j) is the azimuthle angle between p^ 



-jet 



and If J,. 

W+multijet events. In both cases the r coming from W decay has Pi- 



Table 4 summarizes the effect of all these kinematic cuts on the simulated tt and 

1. The R>0.8 cut is 



seen to further reduce these backgrounds to the 25% level. The last two row show the number 
of background events for luminosity £=10 fb~^. Here we have used the LO cross-sections 
for these SM processes for consistency with the SUSY signal. 



Table 5: Events from three SUSY processes for tracker isolation, 



rh 



>3 GeV, 



p^*" >6 GeV. In the last column the total number of events are presented after normal- 
izing the contributions from various sub-processes with their respective cross sections for 



luminosity of lOfb 
R>0.8 cut. 



As in Table 2, the parenthetic entry in each block is obtained with the 



Process^ 


a 
99 


b 
qq + qq* 


c 


a+b+c 
Lumi=10/fb 


No. of event simulated 


300K 


800K 


lOOOK 




^y >250 GeV 

Jet cut 

Me// >750 GeV 

One r-jet 

Mt >50 GeV 


183548 

183548 

155364 

35876 

25083 


624350 
624350 
596914 
114016 
76437 


744353 
744352 
685045 
152093 
105108 




15 GeV< pt <40 GeV 


8065 


25618 


34558 




h 


3209 

(1811) 


17143 
(9681) 


16827 
(9544) 


731 

(414) 


A2 


1991 

(776) 


7043 
(2841) 


8927 
(3651) 


358 
(146) 


xl 


84 
(46) 


10 

(5) 


235 

(112) 


(3) 


xl 


80 

(45) 


57 
(31) 


223 

(130) 


(4) 


w 


2049 
(656) 


873 
(283) 


6084 
(1937) 


191 
(61) 


z 


52 

(22) 


31 

(12) 


230 

(99) 


7 
(3) 


h 


219 
(94) 


396 

(153) 


861 
(343) 


30 

(12) 


r2 


13 

(5) 


1 

(1) 


24 
(14) 


< 1 

(<1) 


None of the above mothers 


368 
(166) 


64 
(30) 


1147 
(508) 


34 

(15) 


Fake as a r jet 
With R cut 


1681 
(40) 


2006 
(62) 


5096 
(126) 


180 

(5) 


Total number of events: Signal 

SUSY Background 

SM Background 








731(414) 
814(249) 
427(110) 
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Figure 1: pri^^ GeV) of softest r-jet for signal(solid) and background processes. 



Table 5 shows the simulated SUSY events for all the three processes of eq.(9). The effects 
of kinematic cuts on the soft r jet events and their distributions to the various sources are 
very similar for the three production processes. The last column gives the total number of 
events from the sum of the three processes normalized by their respective cross sections, for 
a luminosity of 10fb~ . We see from Table 5 that without the R cut the total background 
is almost twice as large as the signal. Imposing the R cut eliminates the fake r background 
and reduces the total background to below the signal size. 

Fig.l shows the pr distributions of the soft r jet signal and the backgrounds before the 
R cut, where the contribution from the fake r jet background is indicated separately by the 
dotted line. The vertical scale on left gives the cross section in fb/0.5 GeV, while that on 
right gives the number of events/0.5 GeV for a luminosity(£) of lOfb^^. Fig. 2 gives the 
corresponding distributions after the R cut of eq. (8). It clearly shows a steep rise of the soft 
r jet signal above the background at the low px end. 

As we see from Tables 4 and 5, the SUSY background is larger than the SM background 
because the former is not suppressed by the selection cuts. Moreover while one can indepen- 
dently estimate the SM background, the SUSY background depends on the relevant SUSY 
masses. It will be very hard therefore to disentangle this background from the SUSY signal. 
It is for this reason that the r polarization plays a very important role in signal extraction. 
The Pr=+1 signal and the Pr = —1 background(SUSY and SM) are known a priori to have 
very distinct R dependences. Likewise one knows the distinct R dependence of the fake r 
background from hadronic jets. Thus for a given pt range of the r-jet events one can make 
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use of the observed R dependence to separate the P^=+l r signal from the Pt=-1 t as well 
as the fake r backgrounds. We shall assume that when real data becomes available, all the 
important SUSY masses, except the small mass difference AM between fi and Xit can be 
roughly estimated (to ~25% accuracy, say), via mass reconstructions using the hard jets 
along with the hard r-jet of pr >40 GeV, as widely discussed in the literature. Then one 
can estimate the size of the SUSY background from this hard r-jet region of the data and 
parametrize it in terms of a cascade decay fit to this data using these SUSY masses. One 
can then extrapolate it to the p^ ""^ <40 GeV region using this parametrization. This will 
provide a fairly reliable estimate of the SUSY background. An observed excess of low p^"''^ 
events over the (SUSY+SM) background, which has been suppressed using the polarization 
cut of R>0.8(Fig.2), will constitute the soft r signal from eq.(3). Thus the procedure for 
estimating the SUSY background from real data is quite clear. In the absence of real data, 
however, we see no better alternative to estimate the SUSY background than the simple one 
adopted here, i.e take it directly from the SUSY event generator along with the signal for the 
illustrative mSUGRA point in the stau-coannihilation region of interest. We feel it suffices 
for our main purpose of demonstrating the importance of the polarization cut in extracting 
the soft r SUSY signal from background. The importance of this cut for signal extraction is 
evident from a comparison of figures 1 and 2. 

The steep slope of the soft r jet signal from the fi decay of eq.(3) over the low pt 
region can be used to extract the signal as well as to measure the tiny mass difference AM, 
responsible for this steep slope. For this purpose, we divide the pr range of the soft r-jet 
of Fig.2 into two parts, px =15-25 GeV and 25-40 GeV. Then we consider the ratio of the 
numbers of events coming from the two parts, i.e 

^' = JT^^. = 2^3 (17) 

A'25-40 = 126 

as a distinctive parameter for extracting the signal and measuring the mass difference AM. 
Being a ratio of cross sections this quantity should be dominated by the statistical error. 
From Fig.2 we can estimate this ratio for the signal+background and for the background 
only. We get 

Ds+B = 467/306 = 1.52 ±0.14; (18) 

Db = 178/180 = 0.98 ±0.14; (19) 

where the statistical errors shown are for lOfb" luminosity run of the LHC. We see that 
this luminosity will be enough to extract the signal at the 3o" level. This will go up to ~10cr 
level at the lOOfb^^ luminosity since statistical error goes down like the v£. 

To estimate the accuracy in the determination of the small mass difference AM between 
fi and Xi, we have chosen a second mSUGRA point near that shown in Table 1, increasing mo 
by 6 GeV. It increases AM from nearly 10 GeV to 15 GeV without practically changing any 
other mass. Fig.3 compares the resulting soft r jet cross sections(signal and background) for 
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Figure 2: Same as Fig.l, but with R cut of eq.(8). 
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The upper sohd and dashed histograms are for signal and background processes for AM =15 
GeV and the lower ones are for AM =10 GeV. These are subject to R cut of eq.(8). 
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the two points, the RHS scale again corresponding to the number of events for a luminosity 
of 10fb~^. One can see a flattening of the low pt slope with increasing mass difference. The 
ratio(18) for the upper curves is 



L)(5+B)(AM= 15 Gel/) = 552/411 = 1.34 ± 0.10. 



(20) 



Thus one can estimate AM to 50% accuracy at ~1.5o' level with a luminosity of 10fb~ . The 
significance level goes up to ha level for a luminosity of 100fb~ , which makes the statistical 
error smaller by factor of 3.2. 

Note that for estimating the small mass differences AM from the pt slope of the soft 
r-jet one has to assume some knowledge of the other SUSY masses. It follows from simple 
kinematics that in the rest frame of the decaying fi the p^""""^ ~ AM. But the boost factor 
relating it to the p^""*^ i'^ ^^^ laboratory frame, pf^/M, depends on the other SUSY masses. 
As mentioned above, we assume that these masses will be known to ~25% accuracy, say, 
from the mass reconstruction programme. Therefore, we need to check the robustness of our 
result to a ~25% variation in the SUSY masses, other than AM. For this purpose we have 
chosen a higher point on the stau-coannihilation strip with tjiq and mi/2 values ~25% higher 
than those of Table 1. The SUSY mass parameter for this point are shown in Table 6. The 
resulting ratio for this point is 



Ds+b{^M = 13 GeV) 



641 

444 



1.44 ±0.11. 



(21) 



Note that this mass difference, AM = 13 GeV, is approximately midway between those of 
the two points represented in Fig.3, AM = 10 GeV and 15 GeV. So we expect the resulting 
ratio (21) to be nearly midway between those of eqs.(18) and (20), which is indeed the 
case. Of course there is room for a more detailed investigation of the sensitivity of the AM 
estimate to the other SUSY masses, which is beyond the scope of this illustrative work. 

Table 6: Masses of superparticles (in GeV), for a higher mSUGRA point in the 
stau-coannihilation region compared to that of Table 1. 



mo 


1^1/2 


tan/? 


9 


Ql 


Qr 


n 


T-2 


xt 


xt 


Al 


X^ 


xt 


xi 


262 


500 


40 


1156 


-1083 


-1043 


219 


425 


390 


635 


206 


390 


622 


635 



It should be mentioned here that we have selected a large value of tan j3 for our illustrative 
mSUGRA point of Table 1 simply because it gives stau-coannihilation region for comparable 
values of mo and mi/2- But we have also checked the results for a low values tan/3 =10. It 
gives stau-coannihilation region for mo << 'n^i/2- However, the size of the soft r-jet SUSY 
signal and background are similar to those of the present analysis. 

Finally let us note that an independent method of probing the stau co-annihilation region 
at the LHC has been investigated in [21] via a ditau-jet signal coming from X2 decay given 
in eq. (11). It contains a hard {pr >40 GeV) r jet along with the soft one, correspodning to 
the r' and r of eq.(ll) respectively. The background is suppressed by taking the difference 
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of opposite-sign and same-sign ditau events. The upper edge of this di-taujet invariant mass 
plot gives an estimate of the tiny mass difference AM, again assuming some knowledge of 
the other SUSY masses. We find in our simulation, however, that the number of such soft 
and hard di-taujet events are suppressed by an order of magnitude compared to the soft 
taujet events analyzed here. Therefore, the present analysis has the benefit of utilizing a 
much larger fraction of the SUSY events, triggered via the selection cut of eq. (14). But, of 
course, the di-taujet channel offers an independent probe of this signal and hence should be 
pursued independently. Let us point out that the r polarization effect can be also exploited 
to improve the result of this di-taujet analysis. In particular imposing the R >0.8 cut of 
eq.(8) on the soft r jet and the complementarity cut on the hard one will ensure that both 
the r jets are leading ones - i.e they carry most of the momentum of the respective r leptons, 
with accompanying soft neutrinos. This will steepen the upper edge of the di-taujet invariant 
mass distribution and improve the resulting estimate of AM, apart from suppressing the soft 
r-jet background. 

5. Conclusion 

The stau co-annihilation region of the mSUGRA model is a region of special interests to the 
SUSY search programme at the LHC In particular one is looking for a distinctive signature, 
which will identify the SUSY signal at the LHC to this region and also enable us to measure 
the tiny mass difference AM between the co-annihilating superparticles, which is predicted 
to be ~5% by the DM relic density constraint. A distinctive feature of this region is that 
a large part of the SUSY cascade decay occurs via fi — > rx?, leading to a r lepton along 
with the canonical missing Exilfrp). Admittedly the resulting r-jets are very soft because of 
the small mass difference between fi and x? states. On the other hand its polarization (P^) 
is predicted to be very close to -|-1. We have shown here with the help of generator level 
Monte Carlo simulation that the positive polarization(Pr=-|-l) of this signal can be exploited 
to extract it from the negatively polarized(PT-=-l) ^-jet background as well as the fake r 
background from hadronic jets. Moreover, the steep pt dependence of the soft r-jet signal is 
shown to provide a distinctive signature for this co-annihilation region as well as a measure 
of the tiny mass difference AM between the co-annihilating superparticles. The significance 
levels of this signal and AM measurement are estimated for the LHC luminosicties of 10fb~^ 
and 100fb-\ 
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